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ABSTRACT 

The next generation of X-ray telescopes have the potential to detect faint quasars 
at very high redshift and probe the early growth of massive black holes (BHs). We 
present modelling of the evolution of the optical and X-ray AGN luminosity function 
at 2 < z < 6 based on a CDM merger-driven model for the triggering of nuclear 
activity combined with a variety of fading laws. We extrapolate the merger-driven 
models to z > 6 for a range of BH growth scenarios. We predict significant numbers 
of sources at z ~ 6 with fluxes just an order of magnitude below the current detection 
limits and thus detectable with XEUS and Constellation-X, relatively independently 
of the fading law chosen. The predicted number of sources at even higher redshift 
depends sensitively on the early growth history of BHs. For passive evolution models 
in which BHs grow constantly at their Eddington limit, detectable BHs may be rare 
beyond z ^ 10 even with Generation-X. However, in the more probable scenario that 
BH growth at z > 6 can be described by passive evolution with a small duty cycle, or 
by our merger driven accretion model, then we predict that XEUS and Generation-X 
will detect significant numbers of black holes out to z ~ 10 and perhaps beyond. 
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1 INTRODUCTION 



The next generation of X-ray satellites promise to be pow- 
, erful probes of quasar activity out to very high redshifts. 
' The European lead project XEU^ is designed to reach a 
sensitivity ~ 100 better than the deepest observations to 
date. Constellation-!X0 would offer at least an order of mag- 
nitude improvement in sensitivity, hopefully paving the way 
for a mission like Generation- which could probe 50 times 
fainter again. 

In this paper we aim to assess the prospects for 
probing the early growth of supermassive black holes 
(SMBHs) at z > 6 with th ese telescopes or simila r fu- 
ture X-ray miss i ons [s ee, e.g., iHaiman fc Loebl |l99^ and 

IWvithe fc Loebl l|2003h for earher work on this]. The ob- 
served correlations between SMBH mass a nd galaxy prop- 
erties , which are now well established (e.g. iGebhardt et al.l 

I2OOOI : iFerrarese fc MerrittI [ioOQj . strongly suggest a tight 
link between the build-up of the stellar mass in galac- 
tic bulges and the mass of the central BH. Most mod- 
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els of the growth history of SMBHs assume that the 
frequent merging of galaxies predicted by CDM-like hi- 
erarchical models of galaxy formati on plays an impor- 
tant role in this proce s s (e.g. Kauffm ann fc Hae hnelt 200^; 



20031 : [ci-oton et all 



Wvithe fc Loeb 2003; Volonteri et al.l 
2006, : Malbon et al. 2006|). Such models, ahhough subject 



to considerable degeneracies, have been reasonably success- 
ful in reproducing the observed luminosity function of AGN 
as well as the inferred BH mass function and their respec- 
tive evolution with redshift. These models generally strug- 
gle, however, to combine the very efficient high redshift 
growth of BHs, suggested by the SMBHs detected at 2 ~ 6, 
with the strong feedback required to reproduce the rapid 
decline i n the fuelling rate of the most massive B Hs at low 
redshift (|Bromlev et al.ll2004l : lMalbon et al.ll2006l ). Further- 
more, at low redshift there is considerable observational ev- 
idence that some, and maybe even most, of the fuelling of 
SMBHs is not connected to major mer gers between galax- 
ies (e.g. iKauffmann et al.ll2003l : iLi et "al . 2007a'l. This may 
be related to the fact that merger-driven models of quasar 
activity have difficulty reproducing the luminosity depen- 
dent density evolution (L DDE) observed for faint AGN at 
redshifts below z ~ 3 (e.g. Uoda et al. 2003: Hasinge r et al.l 
2005: Silverm a n ct a l. 2007; Bongiorno ct al. 2007) [see, e.g., 
Marulli et al.l l|2007h ]. although more sophisticated models 



which follow the simultaneous growth of galaxies and AGN 
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can in corporate feedback in such a way to orchestrate this 
effect (jFontanot et al.ll2006l ). 

We explore predictions for a variety of models for the 
early growth of SMBHs which range from merger driven ac- 
cretion to continuous Eddington-limited accretion. The BH 
growth scenario is crucial in determining the number of ac- 
tive BHs that we expect to detect as faint X-ray sources. 
However, as X-ray observations push to fainter flux limits 
they are expected to discover a mixture of low-luminosity 
objects at low to intermediate redshift, as well as distant 
objects previously at undetectable redshifts. Predictions for 
the redshift distribution of faint sources, which encrypt the 
growth history of SMBHs, therefore require a simultaneous 
study of the faint end of the luminosity function at all red- 
shifts. 

We have adopted a hybrid approach to modelling the 
redshift evolution of the X-ray emission associated with the 
fuelhng of SMBHs. At low redshift {z < 2) we use the ob- 
served X-ray luminosity functions, which are now well estab- 
lished down to very faint luminosities. At intermediate red- 
shifts the faint end of the X-ray luminosity function is still 
subject to some uncertain ty due to the diffi cultly of assigning 
redshifts to faint sources l|Aird et aLlbOOSl l . At 2 < 2: < 6 we 
have therefore adopted a CDM merger-driven model for the 
evolution of the emission from SMBHs, constrained by the 
available optical and X-ray data. The efficiency of BH forma- 
tion, the quasar lifetime and accretion rate, and the quasar 
spectral energy distribution (SED) are all free parameters 
in such a model. The relevant assumptions should obviously 
be guided by our empirical and physical understanding of 
quasars, but there is still much flexibility, particularly near 
the observational limits in redshift and luminosity. 

Observationally, faint AGN are more likely to display 
signature s of obscu ration by gas and / or dus t than brighter 
sourc es (lUeda et al . 2003; Simpsorj I2OO5I : [rrcister et al. 
I2OO6I : iMaiolino et al.ll2007h . Most semi-analytic models of 
quasar activity do not model obscuration effects in suffi- 
cient d eta il to be able to explain these trends. iHopk ins et al.l 
l|2005al lbl. |2006| . 120071 ') have, however, recently presented a 
model for the quasar luminosity function (QLF) at 2: < 5, 
which includes prescriptions for quasar luminosity depen- 
dent fading and absorption, to successfully reconcile the op- 
tical and X-ray QLFs. We follow a very similar approach 
and adopt many of th e same assumptions. However, whereas 
iHopkins et al.l l|2006l ) extract the quasar formation rate from 
the constraints provided by the observed QLF, we use a 
cosmological galaxy mer ger rate from the ex tended Press- 
Schechter formalism fe.g. iLacev fc Coi3ll993l ). This will al- 
low us to extend our model to the redshifts (2 > 6) we 
are predominantly concerned with in a well-motivated way. 
Semi-analytic models are an efficient way of exploring the 
wide range of physically plausible evolutionary scenarios for 
the QLF. As we will discuss in some detail, the faint end of 
the luminosity functions, and therefore the redshift distribu- 
tion of faint X-ray sources, is very sensitive to the assumed 
time dependence of the fading of the emission in different 
wave-bands during galaxy mergers. We utilise the flexibility 
of our model to study the constraints on this fading rate 
provided by the optical, soft and hard X-ray QLFs, and the 
unresolved Cosmic X-ray background (CXRB). 

We then go on to explore several models for BH growth 
at 2 > 6 which are consistent with the observed data at 



< 2 < 6 and span a wide range of possibilities for the rate 
of BH growth at high redshift. 

Our paper is structured as follows. In Section 2 we 
briefly review the constraints on the optical and X-ray QLFs. 
We describe our merger driven model for < 2 < 6 in Sec- 
tion 3, and the high redshift extensions in Section 4.3. In 
the main body of Section 4 we calibrate our model against 
the measured optical and X-ray QLFs. In Section 5 we ex- 
amine the consistency of our models with the CXRB and 
observed logN — logS X-ray source distribution, and present 
our prediction for the source density detectable with future 
missions, before concluding in Section 6. 

Throughout this paper we adopt a cosmological matter 
density Qm = 0.27, baryonic matter density ilb ~ 0.044, cos- 
mological constant Qa ~ 0.73, present day Hubble constant 
Ho = lOO/i km s"^ Mpc"^ = 71km s"^ Mpc"\ a mass vari- 
ance on scales of 8h~^ Mpc as — 0.84 and a scale invariant 
primordial power spectrum (slope n = 1). 



2 OBSERVATIONAL CONSTRAINTS 

We are primarily interested in exploring the constraints on 
the early growth of SMBHs that may be provided by future 
X-ray missions. In order to do this in a conservative way we 
tie our model to a wide range of observational constraints 
on the growth history of SMBHs at low and intermediate 
redshift offered by the QLF in the optical, soft and hard 
X-ray bands, as well as the CXRB. 

The luminosity function of optically selected QSOs is 
well measured up to 2 > 4, with constraints on the bright 
end t hat reach up to 2 ■-^ 6 jCroom et al.l 12004) : I Wolf et al.l 
i; iFan et al.l (20031 . |2004| : Uiang et al.M2007( ). The con- 



I2OO: 

straints on the X-ray QLF extend to 2 ~ 4 at best, but deep 
X-ray surveys surveys reach to considerably fainter limit- 
ing magnitudes than optical surveys. Combined constraints 
of wide-fleld and deep surveys result in a dynamic range 
of ~ 5 orders of magnitude in luminosity at 2 ^ 2 com- 
pared to ~ 3 orders of magnitude at optical wavelengths. 
The l ow luminosities reac h ed by deep X-ray surveys has al- 
lowed IShankar fc Mathuil l|2007j ) to constrain the faint-end 
slope of the optical QLF at 2 ~ 6 from the dearth of X-ray 
detections above 2 ~ 4 — 5. 

The density of optically bright quasars peaks at 2 ~ 3 
with quasar activity decreasing t oward s both lower and 
higher redshift (e . g. Ifiovle et al.1 |200^ : IWoff et al.l l2003l : 
ICroom et aLll2004l : [Richards et al.ll2006l ). The optical QLF 
is well fit by a double power-law with the break luminos- 
ity evolving self- similarly with redshift, although there is 
evidence for the bright end slope flattening toward high 



evidence tor the bright end slope Hattenmg toward high 
redshift JPan et all l20oi : [Richards et al.ll2006l ) [but see also 
IPon tanot et al.l ((2003)]. More recent evidence that the space 
density of lower lumin osities sources peaks at 1 < 2 < 1.5 
(|Bongiorno et al]|2007l ) suggests that the full picture is prob- 
ably more complicated, and that the faint optically selected 
sources show a similar evolution to that seen in X-ray sur- 
veys. The X-ray luminosity function is b est fit with a lumi 
nosity dependent density evolution ( e .g. iMivaii et al.l 



Ueda et al.l l2003l: iFiore et all l2003l : iHasinger et al 



2001 



2005 



Silverman et al.l |2007^ with the source density peaking at 



2 ~ 0.7 — 1. Additional constraints at X-r ay wavelengths 
come from the measurements of the CXRB ( Worslev et al.l 
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l2005l : lHickox fc Markev"itchll2006l ). Obviously the integrated 
emission from faint AGN below the current detection lim- 
its cannot exceed the measured unresolved CXRB and mea- 
surements of the background intensity are thus probably the 
strongest constraints on the faint end of the X-ray luminos- 
ity function. 



3 A HYBRID MODEL FOR THE EVOLUTION 
OF THE EMISSION DUE TO ACCRETION 
ONTO SUPERMASSIVE BLACK HOLES 

3.1 Rationale for a hybrid model 

As motivated in the introduction, we have chosen a hybrid 
approach to modelling the evolution of the emission due to 
accretion onto SMBHs. At low redshift {z < 2) we use the 
well established observed X-ray luminosity functions, which 
extend to impressively faint luminosities. At intermediate 
redshift, where the X-ray luminosity function becomes less 
reliable, we use a merger-driven CDM-like model with a 
range of assumptions for the decline of the accretion rate 
onto the central BH during a merger. This model is cali- 
brated using the extensive optical data as well as the avail- 
able X-ray data. In Section 4 we describe an extension of 
the merger-driven model, as well as models for growth via 
continuous Eddington limited accretion with a range of duty 
cycles , to very high redshift (z > 6). Recentlv lShankar et al.l 
constructed an observationally anchored model for 
the evolution of the SMBH population which indicates that 
the quasar duty cycle increases with increasing redshift, pro- 
viding some empirical impetus for this approach. 

We first discuss the ingredients of our merger-driven 
model in some detail. 



3.2 A CDM-like merger driven model for the 

luminosity function at intermediate redshift 



We assume that 
is predominantly 



the accretion of gas onto SMBIfs 
^ triggered by major galaxy merg- 

ers llKauffmann fc Haehneltl |2000| : iDi Matteo et all l2005l : 
iHopkins et al.ll2005al l. and take the merger rate of dark mat- 
ter haloes in the standard KCDM model for structure for- 
mation as a proxy for the rate of galaxy mergers. This pic- 
ture is admittedly simple but has nevertheless been shown 
to yield a reasonably consistent model for many of the prop- 
erties of the ob served QLF and its evolution at intermediate 
redshifts (e.g. Haohnolt ct al. 1998; Wyithc & Loeb' 2002', 
2003. : Volonteri et al...2003.., ; Croton et al, 2006. : .MaruUi et al., 
20071). 



In a merger-driven model, new quasars continuously 
form at a redshift and luminosity dependent rate. We as- 
sume that the quasars become active with an initial peak 
luminosity Lpcak and then fade according to a "fading law" , 
For a given model for the rate of formation of sources 



dt 
dlogh 



with peak luminosity Lpcak, 'T-(ipcak 
minosity function, ^ 

dn r°° dt ,^ ^ , d^n 



dlogL 



can be written as 



the lu- 



dlogL 



[L, Lpc 



dlogLpcakdt 



dlogLpc3.k- (1) 



We assume that Lpcak is the Eddington luminosity of 
the final mass BH. We further assume that dark matter 



halos host a central BH with mass Afbh given by (e.g. 
IWvithe fc Loebll2003l ) 



Mbh = e(M, z)M, 



q/3 



-1 a/6 



X(l + Z)"/' ( 



(2) 



/2 [ M 



i/3-l 



where Slf; 



n„(i+z)-^ 



1 and 



Ac = IStt^ -I- 82d — 39d^ is the overdensity of a virialised 
halo at redshift z. This relationship is motivated by the ob- 
served correlation between A/ bh and the velocity dispersion 
of th e host galaxie's bulge, a l|Ferraresdl2"oo3 : IShields" et al.l 
120031) . We assume that a may be approximated by Wvir/V2, 
where v^ir is the v irial v elocity of the dark matter halo from 
iBarkana fc Loebl (|200ll ). Empirical estimates of a typically 
fall in the range 4 — 5. We have chosen a = 5 which is con- 
sistent with a simple self-regulated growth scenario in which 
the BH grows until it radiates enough energy to unb ind the 
gas that is feeding it (e.g. ,Silk fc R ees 1998: Haehne lt et al.l 
ll998l : IWvithe fc Loebl [2003h . The merger rate 



dlogLp^^]^dt 

taken as the merger rate of halos with final BH mass such 
that Lpeak = LEdd(Afbh)- Uslug cquatiou ([21 to relate BH 
and host halo mass, we may write 



0.5(M-AM) 



diogLpcakrfi 01 



J AM=0 



=0.25(M-AM) 



N{M, AM, t)dAM, 
(3) 

where the limits of integration ensure that we count only 
major mergers and N{M, AAf, t) is the merger rate of dark 
matter haloes of mass M — AAf and AAf per unit cosmic 
time t, 



N{M, AM, t) 



d^P 



dAMdt 



dn 



d{M - AM) ■ 



(4) 



Here 



is the probability per unit time that a 



halo mass AM will merge with another halo to form a halo 
with mass M from iLacev fc Cold (|l993l ), and ,, ■ . ■ is 
the space density of halos with the appropriate mass differ- 
ence from Press fc Schechterl (Il974 ) (with modifi c ation from 
gheth fc Tormcn' ('1999'). see iRhook fc Haehneltl (j2006l ) and 
Wyithc fc Locb (2003) for similar calculations). 

Equation ((3} becomes inaccurate as -j^^ approaches 
the Hubble time at the relevant redshift. This could in 
principle be relevant at high redshift, however at high red- 
shift we are mainly concerned with the bright end of the 
QLF for which the fading time-scale remains short (see Sec- 
tion |3]3l). Note that for ^^^^^ (L, -Lpcak) ~ constant (corre- 
sponding to an exponential light-curve) or ^ (L, Lpeak) oc 
5{L ~ LLpcak) (corresponding to a top-hat light-curve) the 
shape of the QLF is identical to that of n(Lpcak). 



3.3 Gas accretion, quasar lifetime and fading 
laws, and the minimum black hole mass 

As discussed in the previous section, the faint end of the 
QLF is very sensitive to the assumed fading law. When 
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two gas-rich galaxies merge, the resulting tidal torques drive 
large amounts of gas into the central region, providing fuel 
for the rapid growth of a SMBH. The fading of the quasar 
emission due the gas accretion onto a central SMBH is gov- 
erned by the infall of gas to the centre and its subsequent ac- 
cretion onto the SMB H during the late stages of the merger. 
iHopkins et all l|2005ah used numerical simulations of galaxy 
mergers, with a prescription for t he subsequent a c cretio n 
of gas onto a central SMBH from iDi Matteo et all (l2005f). 
to obta in a physically motivated fading law. iHopk ins et al.l 
l|2005al ) present a power-law representation for this fading 
law, , , which has been fit to the results of several hun- 

' dlogL ' 

dred simulations of mergers between equal mass galaxies 
resulting in quasars with peak bolometric luminosities be- 
tween 10* and 10^^ Lq. The power-law index depends on 
Lpeak, with higher luminosity quasars expected to spend rel- 
atively more time close to their Eddington limit. 



dt{L, Lpeak) 

dlogL 

(Lpeak) 



109 L, 



= -0.95-^0.321ogio(Lpeak/lo'' Le 



(5) 
(6) 



where tg « 10^ yrs and ol > -0.2. iHopkins et"al] (j2005al l 
model the quasar luminosity over the entire merger pro- 
cess and therefore equation ((5} includes the activity of the 
quasar prior to the luminosity reaching Lpcak- This distinc- 
tion should not be important for our modelling. 

In order to investigate the effect of the quasar light- 
curve on the the faint end of the lumi nosity function, we con - 
sider the fading law as proposed by IHopkins et ahl (|2005al ) 
as well as a fading law where quasars fade more rapidly and 
spend equal amounts of time at each luminosity (relative to 
the peak luminosity) independently of Lpeak, 



dt 



dlogL 



L 



peak 



(7) 



Here q_l is a non-positive constant and larger values of |aL 
result in light curves for which quasars spend relatively more 
time at sub-Eddington luminosities. Note that values of q:l 
that are much less than zero are excluded by the observa- 
tional data. 

The faint end of the luminosity function is also sensitive 
to the minimum mass of a BH powering a quasar, Mbh,min, 
which for Eddington limited growth is equivalent to the min- 
imum peak luminosity of a quasar, Lpeak, min- There are very 
few claimed detec tions of SMBHs with masse s smaller than 
lO'^Mp) (but see iGreene fc Hoi l2007al ). and iGreene fc Hoi 
l|2007bl ') find that the mass function of local low-luminosity 
active BHs turns over at Mbh ~ lO'' M0. It is also a mat- 
ter of intense debate whether the build-up of SMBHs by 
hierarchical merging extends to intermediate mass BHs in 
the mass range 100 M© - 10^ Mq (e.g. iHaehnclt 2004; 
IVolonteri et~al]|2007l ). We therefore also explore the effect 
of varying the minimum BH mass, Mbh, min, in our models. 



3.4 Modelling the Spectral Energy Distribution 

To relate the intrinsic bolometric QLF to observational data, 
we also require a model for the intrinsic SED and a pre- 



scription for the relative obscuration of quasars at different 
wavelengths. 

The analyses of samples of optically selected (type-I) 
quasars suggest that there is a luminosity dependent corre- 
lation between the X-ray and optical luminosity, which may 
be described in terms of the power-law index connecting 
the optical and X- ray flux, aox ■ ceox may be written (e.g. 
ISteffen et aLll2006l ) 



aox = 



logio(L^(2500 A)/L^(2 keV)) 
logi(,(;/(2500 A)/u{2 keV)) ' 



/ m I ' (8) 
y ergs/s/Hz / 

where A,B > corre sponding to a n ant i-correlation of 
ctox with optical flux. 'Ste ffen et al.l (|2006| ) estimate A — 
0.137 ± 0.008 and B = 2.638 ± 0.24 for an optically selected 
sample of qua sars with no appare nt intrinsic absorption or 
radio activity. ISteffen et al.l (|2006l ) find no strong evidence 
for evolution of this relationship with redshift, however the 
data is consistent with aox decreasing slowly with redshift 
as ~ 0.012. 

We a dopt a modified version of the model SED pre- 
sented in IHopkins. Richards, fc Hernguistl (|2007l ) to calcu- 
late the (luminosity dependent) SED for a given observed 
frequency, i^obs, quasar redshift, z, and bolometric luminos- 
ity, Lboi. This model assumes the SED to be essentially 
a broken power-law which mimics reprocessed emission in 
the NIR, the UV bump and the hard X-ray excess features 
with an exponential cut-off at 500 keV. The SED is nor- 
malised so that the observed relation between optical and 
X-ray emission described by equation ([8)l with A = 0.109 
and B — 1.739 is preserved. 

Note that IHopkins et all l|2007h found that the lumi- 
nosity dependence of the relation between optical and X-ray 
emission is more important for reconciling the optical and 
X-ray QLFs at < z < 5 than the detailed modelling of 
the features in the SED. The anti-correlation of aox with 
Ly(2500A) results in a narrower spread of X-ray luminosi- 
ties compared to the corresponding optical range. The model 
X-ray QLF is thus steeper than the equivalent bolometric 
QLF. 

We find that the bright end of our model X-ray lumi- 
nosity function becomes too steep compared to the measure- 
me nts in both the soft a nd hard X-ray bands when we use 
the I Hopkins et al] |2003) SE D. Including a dispers ion in the 
bolometric corrections, which IHopkins et al.l (|2007| ) note has 
the largest effect on the bright end of the X-ray QLF, would 
smooth the QLF and somewhat alleviate this. We have cho- 
sen instead to alter the luminosity dependence of the SED 
model to fit the data as described below. 

To achieve a good match with the bright end of the ob- 
served X-ray QLF we adopt an SED for Lboi > 10^* L© with 
A = 0.104 and B = 1.739. We use the same power-law in- 
dices, break-points and exponential cut-off as IHopkins et al.] 
(|2007h but have not modelled the reflection component or 
the NIR emission. This choice of parameters increases the 
power in the soft and hard X-ray bands by approximately 
40 per cent at Lboi ~ 10^* L©, increasing to 45 per cent 
at 10^^ Lq. The power at UV wavelengths changes mini- 
mally. To ensure a s mooth QLF we take an average of the 
iHopkins et al] l|2007h and the above modification for bolo- 
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logio[(t-to)/years] log , o[(t-to)/yea rs] 

Figure 1. Light curves used to generate model QLFs for quasars with intrinsic, peak, bolomctric luminosities of 10^^ (solid line), lO'^^ 
(dashed line) and lO"*^^ (dot-dashed line) Lq. The left panel shows the slow fading law, corresponding to equation JSj with the average 

lifetime (tg) increased by a factor of two. The right panel is for the rapid fading law ^I'^^j^ = 1-7 X 10 yrs (77-^ — j corresponding 
to rapid (almost exponential) luminosity independent fading. 



metric luminosities between 10^^ and 10^* Lj 
the logarithmic bolometric luminosity. 



weighted by 



3.5 The effect of absorption on the SED 

Light emitted from the accretion disk of the AGN is repro- 
cessed during transmission through its host galaxy, the IGM 
and our own galaxy. At the X-ray energies and redshifts we 
are mainly interested in (observed frame 0.5 — 10 keV) the 
absorption of X-rays by neutral hydrogen in the IGM can 
be neglected. Galactic absorption is generally corrected for, 
and therefore we are mainly concerned with absorption at 
the quasar redshift. Assuming that emission at all wave- 
lengths is obscured by the same body of gas, the absorption 
at a particular wavelength can be determined given a column 
density distribution of neutral hydrogen (which may depend 
on parameters of the system) a dust-to-gas ratio and a red- 
dening law for the dust absorption. W e again adopt similar 
assumptions as iHopkins et all l|2005bl ) for our model, which 
we discuss in the following two sections. 



exp{-[logio(7VH/iVH)]V(2<7^H)} 



0"JV„ 



10 cm 



1.2. 



(9) 



Note that the mean column density, iVn, and dispersion, 
(Tnh, are predominantly determined by the instantaneous 
bol ometric luminos it y of th e quasar L. 

liO |2005al ) also demonstrated that the 



iHopki 



ms et 

shape of this distribution for quasars with a B-band luminos- 
ity above lO'^^ Lq is consistent with a statist ical analysis of 
the reddening found toward SDSS quasars in lHopkins et al.l 
(2004). This analytic fit is derived from simulations of 
gaseous disc galaxies, and of course may not be represen- 
tative of the typical accretion rate onto black holes at all 
redshifts. However at the redshifts that we are applying it 
(z > 2) galaxies are expected to be gas rich. 



3.5.1 Column density distribution of the absorbing gas 

The column density distribution of absorbing gas is ex- 
pected to be related to the fading law; gas both feeds the 
BH and obscures the radiation released. However the as- 
trophysical relationship is complicated and depends on, for 
example, the surroun ding star formation an d the geometry 
of the accretion disc. iHopkins et all l|2005al ) have used the 
numerical simula t ions o f quasars fuelled by gas accretion by 
iDi Matteo et al.1 (|2005l ) to construct a parametrised lumi- 
nosity and time dependent model for the obscuring column 
density. The probability P that a quasar is obscured by a 
given total hydrogen column density A^h depends on the 
amount of time a quasar remains in a given accretion phase. 
The simulated data is fit with a log-normal distribution 



3.5.2 Dust and gas absorption 

Since we are predominantly comparing our model to optical 
data quoted in UV magnitudes we hav e used the redden- 
ing law presented in lGaskell et al. I (|2004l ) to model dust and 
gas absorption. The reddening law is normalised such that 
the optical depth in the V-band (5500 A) i s the sam e as 
for an SMC-like reddening curve [taken from (|Peilll993 )] for 
a galaxy with the metallicity of the Milky Way. This re- 
sults in a V-band opti cal depth identical to that adopted by 
iHopk ins et al. I (|2005bl ) who assumed a SMC-like reddening 
curve with a gas-to-dust ratio as for the Milky Way. At X-ray 
wavelengths the absorption is dominated by photo-electric 
absorption and Compt on scattering by hydrogen. Following 
iHopk ms et al.l (|2005bh we assume an average neutral gas 
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fraction of 0.35 and assume that the ionised component con- 
tributes to the photoelectric absorption, but not the optical 
reddening or Compton scattering. We note that this model 
makes the simplistic assumption that the distribution of gas 
and dust is spatially uniform. 



3.6 The model luminosity function taking 
absorption into account 

We now turn to calculating the observed luminosity function 
for a given redshift z and observed wavelength v (or more 
usefully for comparison with X-ray observations, wave-band 

Given the probability that a quasar with a given in- 
trinsic bolometric luminosity L' is obscured by a given col- 
umn density, the probability that the intrinsic specific lu- 
minosity L'^ is observed to have a specific luminosity L^, 
P(L„\L'^)d\ogLl , may be written 



P{L,\Ll)dlogLi 



-log 



= P{Nh = —log 



■rflogl/L, 



dlogL'„ 



where is the absorption cross section at u in units cm'^. 
The absorbed luminosity function at an observed wavelength 

!/ and redshift z, ^i^^j^ \ {z), can be calculated by inte- 



grating over the unabsorbed bolometric QLF 



dti 



dn 



dlogLi, 



{z) = 



dn 



iog{L' )=iog{L„(i+,)) dlogL^^j^^^^ 



(10) 



The convolution is identical if we want to consider the 
luminosity function in band (uaji^b), except the change of 
variable in equation (|10|l becomes '^'"^i"-'^-" 



Here < cr^ > is the average absorption cross section in the 
band, weighted by the observed luminosity vL^e'""^" , 



(11) 



3.7 The blowout phase 

The probability distribution of the absorbing hydrogen col- 
umn density given in equation ([5} will not be an adequate 
description when the luminosity is close to the peak lumi- 
nosity, as much of the gas is expecte d to be blown a way by 
radiation pressure. As discussed by iHopkins et al.l (|2005al ) 
the true column density distribution in the simulations is bi- 
modal, with quasars in the last e-folding of growth having a 
much lower obscuring column density distribution. This so 
called "blow-out" phase lasts approximately 10 per cent of 
the tot al time the quasar spends accreting. iHopkins et al.l 
lj2005al ) suggest that this bi-modality can explain why opti- 
cally selected quasars are observed to ha ve lower obscurin g 
column densities for brighter quasars (e.g. lUeda et al.ll2003h . 



counter to the positive correlation between L and Nh in 
equation (|9}. 

We have taken a simple approach to modelling such a 
blowout phase; we assume that a fraction /blowout = 0.1 of 
quasars within a factor of e of their peak luminosity expe- 
rience no intrinsic absorption. The resulting absorbed lumi- 
nosity function is then larger than that in equation (|10p by a 
term -^{^Jl^ > i ) x /blowout , where -j^{^J^l) is 



the QLF of sources with specific luminosity that have to- 
tal bolometric luminosity within e of their peak luminosity. 
diogL ("l — ^ )i increases with Lboi, and approaches the 

original (without blow out phase) value of ^i^gj^ for large 
Lboi. This prescription for the blowout phase essentially puts 
a lower limit on the ratio of the absorbed to unabsorbed lu- 
minosity function at high luminosities equal to /blowout. 

The parameters of the blow-out phase affect the rela- 
tive number of bright optical and X-ray luminous sources. 
In particular, if we didn't include a blow-out phase (and 
were therefore assuming that more of the bright optical 
sources are obscured) then we would predict more X-ray 
bright sources. However we find that the absorption in the 
soft X-ray band behaves similarly to that in the optical and 
therefore if we did not include the blow-out phase it would 
be very difficult to reproduce the space density of quasars 
that are bright in soft X-rays. 



4 EVOLUTION OF THE MODEL 

LUMINOSITY FUNCTION IN DIFFERENT 
WAVEBANDS AND COMPARISON WITH 
OBSERVATIONAL DATA 

4.1 Calibrating the merger-driven model for the 
evolution of the luminosity function at 
intermediate redshift with observations 

We vary the parameters {tq, eo, Lpcak.min) in order to obtain 
an acceptable fit to the observed luminosity functions at 
rest-frame 1450 A, and in the soft and hard X-ray bands. 
We determine the normalisation of our models (governed 
by tq and eo) by comparing to the observed optical QLF, as 
these constraints span the broadest redshift range. However, 
as we will discuss below, the optical data has little power to 
constrain Lpcak.min and constraints on this value come from 
low-luminosity X-ray observations alone. 

We find that assuming the luminosity independent fad- 
ing law in equation (O, it is difficult to accommodate val- 
ues of az, much less than zero, due to the intrinsic steep- 
ness of the cosmological merger rate we are adopting for 
the quasar formation rate. However for ql = —0.01, tq — 
1.74 X 10*^ years and eo ~ 10"^"^ we recover a good fit to 
the optical QLF at z = 2 — 6 (see Figure 2). This value of a_t 
corresponds to a light-curve for which the luminosity drops 
off almost exponentially with time, we will therefore refer to 
this model as the "rapid fadin g" model. The correspo nding 
light curve is compared to the iHopkins et all (|2005al ) fight 
cur ve (with some a djustment described below) in Figure 1. 
The lHopk ins et al] (2005a) fading law, with appropriate eo, 
also provides a reasonable fit to the optical data at z = 2 — 6 
as described by equation ^ . To put it on equal footing with 
what we have done for the fit with the rapid fading law 
we allow the characteristic time, tg, to vary. We find that 
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Figure 2. Model rest-frame 1450 A QLF at z = 2,4,6,8,10,12 for a merger driven model with a slow fading law {solid line) and the 
rapid fading law (dashed line) as described in the text. Note that we have neglected absorption by neutral gas in the IGM which will 
affect what is observed in practice at z > 6. The hollow circles are the measurement of the faint end QLF from the CO MBO-17 survey 
(Wolf et al. 2003). The black triangles at 2 = 2,4 show the luminosity functi on measured bv [Bqngiorno et al. 1 ll2007t) in the redshift 
interval 2.1 < z < 3.6. The small black diamonds at 2 = 2 are the results from lCroom et aL ] ||2004| ') using the SED from [Hopkins et al.l 
to convert from B-band to rest frame 1450 A luminosities. The grey bow-ties represent t he constraints on the bright end QLF 
measured in the SDSS survey (Fan et al. 2002, 2004). The thick black line at z = 6 is the fit from I Jiang et al.l ||2007| ). The thin grey line 
at z = 6 is the Shankur & Mathur (2006) constraint on the slope of the QLF from the dearth of X-ray sources. 



the fit is improved when tg is increased by a factor of 2 to 
2 X 10® years. A shghtly larger value of eo = 10~*'®^ than 
for the rapid fading model is required to offset the smaller 
amount of time that bright quasars spend at their Eddington 
luminosity. We refer to the model with this fading law as the 
"slow fading" model. Since the two models provide compa- 
rable fits to the optical data, the optical data alone appears 
to offer little power to constrain the luminosity dependence 
of the fading law. 

Note that there is an excess in the predicted number of 
optically bright quasars at 2 = 2 for each fading law. The 
most plausible explanation for this discrepancy between our 
model and the data is probably our neglect of AGN feed- 
back. A GN inject large amounts of heat into the ir surround- 
ings (e.g. lDunn fc Fabianll2006l : lBest et al.ll200iD and should 
be capable of suppressing cooling flows in dark matter ha- 
los with masses above ~ few x 10^^ Mq - resulting in 
the formation of groups and clusters of galaxies rather than 
sup er-sized quasars in very large dark matte r halos (see, 
e.g. ISiiacki fc Springel|[2006l : ISiiacki et al. I l2007l : iRines et all 
l2007i r This effect is expected to be more prominent at low 



redshifts when the typical masses of merging dark matter ha- 
los is largest. Indeed most semi-analytic models of the quasar 
population require some arbitrary h igh mass cut-off in order 
to fit to the low redshift da ta (e.g. iKauffmann fc Haehneld 
I2OOOI : ICattaneo et al.ll2007l) . To capture this effect in our 
model in a redshift dependent way would be difficult and 
we do not attempt this here. It is also for this reason that 
our model, and similar semi-analytic models for the QLF, 
do not reproduce quasar luminosity function below 2 ~ 2 
very well. 

With tq and eo chosen to reproduce the optical lumi- 
nosity functions we find fairly good agreement with the ob- 
served 2 ~ 2 QLF in the hard X-ray band, and consistency 
with the constraints on the 2 ~ 4 hard X-ray QLF (see Fig- 
ure 3). The value of the minimum BH mass (or minimum 
peak luminosity) does not affect the QLF at the luminosi- 
ties relevant for optical constraints, however the constraints 
on the X-ray QLF are approximately 2 orders of magnitude 
deeper at 2 = 2 than in the optical. We find that a cut-off 
at Lpoak = 10^^ Lq is required for the rapid fading model 
in order to avoid over-predicting the number of sources in 
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Figure 3. Model rest-frame 2 — 8 keV QLF at z = 2, 4, 6, 8, 10, 12 for a mer ger driven model wi th slow fading [solid line) and rapid 
fading (dashed line) as in Fi gure 2. The filled g rey circles are estimates from lBarger et al and the hollow grey ci rcles the rest 

frame 2 — 10 keV data from lUeda et al. I 1 120031) . The thick grey line in the 2 = 4 panel corresponds the estimate from iBarger et al.l 
(I2OO5I') assuming that all spectroscopically unconfirmed sources lie in the z = 3 — 5 bin. The black markers are the hard X-ray QLF from 
[Silverman et al ] iHoOTI); the hollow diamonds for 1.5 < z < 2, the hollow triangles for 2 < z < 3, the hollow squares for 3 < 2 < 4 and 
the fil led circles for 4 < 2 < 5.5. The stars with horizontal and vertical error bars are the measurement for 2.5 < 2 < 3.5 from lAird et al.l 
1I2OO8I '). 



the lowest luminosity bin in the hard X-ray QLF at z = 2. 
We apply the same cut-off in Lpeak to the Hopkins-fading 
model and find it over-predicts the rather shallow faint end 
of the observed luminosity function at z = 2 in the hard X- 
ray band (see Figure 3). Note that an Eddington luminosity 
of 10^^ Lq corresponds to a BH mass ~ 3.10^ Mq and there 
is certainly evidence f or the existence of B Hs with masses 
smaller than this (e.g. iGreene fc Hdl2007al ). 

In Figure 4 we plot our model observed fra me 0.5—2 keV 
QLF and compare it to the results of Miyaji et al.l (|2000l ) , 
iMivaii et all (120011 ) and iHasinger et all l|2005l ). Even with 



the increased power at X-ray wavelengths for the brightest 
sources, we still seem to under-predict the space density in 
the brightest bins at 2 = 2 and 4. Note, however that the 
highest L data point is at each redshift is calculated from 
only once source with an extremely high soft X-ray luminos- 
ity of ~ 10^* Lq, and therefore this discrepancy is probably 
not significant. 

We find that whilst the normalisation determined from 
the optical constraints reproduces the counts well at in- 
termediate luminosities, our models consistently overpre- 
dict the number of faint sources at z = 4, and to a lesser 



extent at 2 = 2, w hen compared t o the type-1 AGN lu- 
minosity function of IHasinger et al.l (|2005l ) . This inconsis- 
tency may be partiall y attributable to the restriction of the 
IHasinger et al.l (|2005l ) to type-1 (unobscured) AGN and in- 
completeness of the optical identification. If we take the 
data at face value, then given that the merger-rate of galax- 
ies steepens significa ntly toward higher redshifts, then the 
IHasinger eraLr(l2005l ) data suggests that for a merger-driven 
model of quasar activity Z/pcak,min must be increased by ap- 
proximately an order of magnitude to 10^^ Lq . Such a cut-off 
brings the rapid fading model into agreement with the hard 
and soft X-ray data at z = 2, 4, and eases, though does not 
eradicate, the excess of faint objects predicted by the slow 
fading model. This modification would not alter the agree- 
ment of either model with the optical constraints. 

Invoking a cut-off at Lpeak.min ~ Lq corresponds 
to a minimum mass of active black holes of ~ 3 x 10^ Mq 
in our model. Whilst the mass fun ction of local BHs m ay 
turn over at around this value (e.g. I Greene & Ho 200J7a||B). 
this trend is highly uncertain and direct measurements of 
BH masses at higher redshift are not yet available. 

However, it appears possible that the current surveys 
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Figure 4. Observed frame 0.5 — 2 keV QLF at z = 2, 4, 6, 8, 10, 12 for the same models as in Figure 3. The thiclc grey line is the best-fit 
luminosity dependent density evolution m odel from (|Mivaii et al.l2o"ooh (LDDE2) . The hollow mark ers are the binned lum inosity function 
for z = 1.6- 2.3 from lMivaii et al.l 1I2OOII '). The black triangles denote the type-I AGN QLF from llHasineer et al.ll2005^ . 



have u nderestimated th e space density of faint sources at 
z = 4 (jAird et aLll2008l ). Note that if the measurements of 
a small space density of faint hard and soft X-ray selected 
quasars a,t z — 4 and the sudden drop in the density of 
X-ray quasars between z = 2 and 4 consolidate, then the 
density of high redshift X-ray quasars, and therefore the 
number of sources detectable by the next generation of X- 
ray satellites, may be lower than we are predicting with our 
fiducial merger-driven model. In the discussion that follows 
we will also discuss results for a slow and rapid fading models 

with Lpcak.min = 10^ Lq. 

4.2 Extrapolating the evolution of the luminosity 
function in the merger-driven model to very 
high redshift 

A prominent feature in our models of the high redshift hard 
and soft X-ray QLFs is the break at low luminosities in the 
rapid fading model due to the cut-off Lpoak.min = lO'^^ Lq. 
No such break is observed for the slow fading model where 
the sources with Lpoak ~ ipoak.min fade gradually. Above 
the break, the X-ray QLF exhibits little dependence on the 
fading law. This is because at high redshift the dependence 
of the merger rate on luminosity in our model results in a 
luminosity function with steeper slope than the fading rate 
of sources for either fading model. There exists therefore a 



degeneracy in the prediction of the number of faint sources 
between the slope of the assumed fading law and the mini- 
mum peak luminosity. Values of Lpoak, min < 10^^ L0 are cer- 
tainly plausible at high redshift; theoretical models for the 
formation of the first BHs suggest that seed black holes with 
mass es between 100 and 10"* Mp) may form as early as 2; = 20 
(e.g. iMadau fc Ree3l200ll : IVolonterillioO^ . In Section 6 we 
explore the possibility of a lower value of -Lpcak,min = 10^ L© 

(Mbh.min ~ 3 X 10-* Mq). 

4.3 Alternate models for the early growth of 
supermassive black holes 

The mechanism for the growth of very high redshift BHs 
is very uncertain. Estimates for the formation of the first, 
or seed, BHs vary widely. BHs with masses 10^ — 10^ Mq 
may form from populatio n-Ill remnants at z = 20 or ear- 
lier (|Madau fc Reesll200ll ). To grow BHs as large as a few 
X 10^ Mq by z ~ 6 from a pop-Ill remnant mass seed 
requires more or less conti nuous accretion if the accretion 
rate Eddington limited (e.g. lArchibald et al. l l2002t ). lLi etHI 
(20 07b ) find that in a merger-driven scenario, BHs large 
enough to power the brightest 2 ~ 6 quasars may be grown 
from seed BHs with mass 10^ M©. In this model most pro- 
genitors in the (hierarchical) merging history of the SMBH 
grow at sub-Eddington rates for the majority of the time. 



10 K.J. Rhook & M.G. Haehnelt 



The overall growth, however, is dominated by the fastest 
growing BHs at each step in the hierarchy. In aggregate the 
growth is similar to almost continuous Edd ington limited 
accret ion of a single seed BH. Following iRhook fc Haehneltl 
l|2006l ) we consider a passive evolution scenario in which a 
fixed comoving density of BHs evolve by accreting at their 
Eddington limit with a fixed duty cycle /duty In this model 
for the early BH growth, the ensemble average luminosity of 
a BH increases with decreasing redshift as, 



L 



L{z = 6)e 

C(J e ^acc 



5 X 10* yrs—^ 



(12) 
(13) 



where (Je is the Thompson scattering cross section. 

The exponential growth of BHs in this model makes it 
unsuitable for modelling the QLF below 2; ~ 6. We therefore 
consider this model as a possible alternative to the merger- 
driven evolution only at z > 6. We adopt the rapid fading 
merger-driven model with Lpeak.min = 10^ Lq at z < 6. 

Extending the absorption laws used in the previous sec- 
tion, which were determined from estimates of the average 
column density along the line of sight to a pair of merg- 
ing galaxies, is probably not meaningful for continuously 
accreting BHs. We therefore simply assume that a fixed, lu- 
minosity independent, fraction of the accreting BHs are un- 
obscured, and that the rest are Compton-thick (completely 
obscured in the optical and X-ray) . The total duty fraction 
is then made up of the duty fraction of BHs that are ac- 
creting in an X-ray luminous phase fium and those that are 
accreting in a Compton thick phase /obsc. 

Note that for a model in which the BHs are growing 
continuously (although with some duty cycle) the concept 
of a fading law is also no longer meaningful since we are 
assuming that there is always enough gas for the BHs to 
accrete at their Eddington limit. 

We have explored the evolution above z = 6 for two val- 
ues of the total duty fraction, /duty = 0.1, 1.0. The fraction 
of obscured sources at high redshift is unkno wn but is cer- 
tainly significant at lower red shift (e.g. .Risaliti et al.lll999l : 
iMartfnez-Sansigre et al.ll2007l ). We assume that half of the 
sources are totally obscured in each case, which is consis- 
tent with the constraints from population synthesis models 
of the CXRB (|GiUi et al.1120071 ). 

For our models with Eddington limited BH growth at 
z > 6 our assumptions governing the accretion rate and in- 
trinsic absorption change discontinuously at 2 = 6 in such a 
way that the observed hard X-ray QLF changes smoothly. 
This alternate model for the early growth of SMBHs should 
thus be regarded only as a demonstration of the uncertainty 
in the growth rate of BHs, and therefore the number of de- 
tectable X-ray quasars, at very high redshifts. 

In Section fS.SI we explore the predictions for the number 
of quasars detectable at X-ray wavelengths for both passive 
and merger-driven BH growth scenarios. We first compare 
our merger-driven models to the constraints from the ob- 
served fiux distribution of X-ray sources and the CXRB. 



5 X-RAY NUMBER COUNTS AND REDSHIFT 
DISTRIBUTION AT FAINT FLUX LEVELS 
AND THE EARLY GROWTH OF 
SUPERMASSIVE BLACK HOLES 

5.1 X-ray source counts and redshift distributions 

The main aim of our modelling is to predict the X-ray source 
counts and the corresponding redshift distribution at faint 
fiux levels. 

The X-ray source counts ("the logA*' — \ogS relation- 
ship") are well described by a double power-law, with the 
slope fiattening for fluxes below ~ 1.5 x lO'^* erg s~^ in 
the soft X-ray band and below ~ 4.5 x 10^^ erg s~^ in 
the h ard X-ray band (e.g. ICowie et al.ll200^ : iMoretti et al.l 
l2003h . More recently, the constraints on the number counts 
have been extended to lower sensitivities by fluctuation anal- 
yses of the unresolved background. These constraints gen- 
erally still allow for an upturn in the logA^ — logS* rela- 
tionship below the detectio n limit for resolved sources (e.g. 
iHickox fc Markevitchll2006l ). 

In the top panel of Figure 5 we plot the source counts 
as a function of limiting fiux sensitivity in the hard X-ray 
(2 - 10 keV) band and in Figure 6 for the soft (0.5 - 2 keV) 
band. The left panels show the results for the slow fad- 
ing model, the middle for the rapid fading model with 
ipoak.min ~ 10^^ Lq and the right for the rapid fading 
with Z/poak,min = 10^ Lq. Iu each case th e over-plotted data 
is tak en from the fiuctuation analysis of lMivaii fc Griffithsl 
(| 2002| ) for the Chandra Deep Field North (CDF-N). Recall 
that below z = 2 we have supplemented our model for the 
X-ray QLF with a fit to the observed QLFs (calculated at 
discrete reds hifts z = 0.01 , 0.6, 1 .2). We have used the PLE 
model from iBarger et al.l (|2005|) for the hard X-ray band 
and the LDDE2 model from l|Mivaii et aLllioOlh for the soft 
X-ray band. Therefore by construction our models should 
reproduce the observed source counts if our model is a rea- 
sonable approximation at z ^ 2. 

Our models are consistent with the constraints on the 
hard X-ray counts for both the rapid and slow fading laws. 
However, for the soft X-ray counts only the rapid fading 
models are consistent with the observed counts. Using the 
slow fading law our model slightly overshoots the faint soft 
X-ray counts and lies above the limi t sugg est by the fluc- 
tuation analysis of lMivaii fc Griffit hs' (2002*) , as is expected 
from the inconsistency with the measured QLF at z ~ 2 
discussed in Section 4.1. For the slow and rapid fading mod- 
els with Lpoak.min = lO^'^ Lq , the uumbcr counts below the 
current detection levels are reduced, bringing the slow fad- 
ing model into agreement with the fiuctuation analysis and 
resulting in a fiat source distribution at low fiux levels for 
the rapid fading model. 

In the lower panel of Figure 5 we plot the normalised 
redshift distribution for the source counts above three 
fiux/sensitivity levels. The grey line is for the approximate 
sensitivity of current X-ray satellites (3 x 10~^^ ergs/s/cm^), 
the solid black line for an order of magnitude fainter (3 x 
10~^^ ergs/s/cm'^) and the dotted line for the goal sensitiv- 
ity of a 1 Ms observation with XEUS (3 x 10"^** ergs/s/cm^, 
but see the discussion in Section r5.3|) . We note the strong de- 
pendence of the redshift distribution on the fading law. For 
the slow fading model the probability for source detection 
peaks at 2 ~ 1 — 3 for each of the chosen flux limits, with 
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Figure 5. Results for the expected number counts of sources in the 2 — 10 keV band for the model with slow fading (left) and rapid fading 

with i/peakjinin 

= 10" L0 (middle) and 

^peak,min — 10^ ^^0 (lis^^^)- Note that the models have been supplemented with observational 
QLFs below z = 2 as described in the text. The top panel shows the counts integrated from z = 0,2, 4, 6, 8, 10 to 2 = 20 as labelled. 
The hollow markers with error bars (bow-tie) are the estimate of the resolved (unresolved) source density in the Chandra deep field by 
Miyaji & Griffiths (2002). The bottom panel shows the redshift distribution of the counts for sources above flux levels of 3.e — 18 {dotted 
line), 3.e — 17 {solid line) and 3.e — 16 {grey line) in cgs units. 




Figure 6. Results for the expected number counts in the 0.5 — 2 keV band for the models in Figure 5. At z < 2 the models have been 
supplemented with the LDDE2 fit to the measured 0.5 — 2 keV QLF from Miyaji, Hasinger & Schmidt (2000). The hollow markers with 
error bars (bow-tie) are the estimate of the resolved (unresolved) source density in the Chandra deep field by Miyaji & Griffiths (2002). 



12 K.J. Rhook & M.G. Haehnelt 



the peak shifting to slightly higher redshift for the lower flux 
limits. For this model the source counts remain dominated 
by fading objects at intermediate redshifts where the space 
density of quasars is at its peak. 

For the rapid fading models the source distribution is 
sharply peaked at z ~ 1 for the brightest flux limit. This is 
not surprising: this is where most of the known sources lie 
and we have used the observed data at low redshifts as our 
model QLF. At fainter flux limits a second broader peak 
around 2 ~ 4 emerges for the rapid fading models. The 
broad high redshift peak at faint flux level in our model is 
due to the steepening of the merger rate at z > 2 combined 
with the increase in the comoving volume element out to 
z ~ 3.5. 

The emergence of two distinct populations at faint flux 
levels, at low and high redshift, can be explained physically 
by an evolution of the Eddington ratio and the characteris- 
tic lifetime with redshift. Observationally, optically selected 
high redshift and/or high luminosi ty AGN display near Ed- 
dington accretion rat es (see, e.g. iMcLure k, Dunlopl 12004 
iKoUmeier et al.ll2006l ). whereas hard X-ray selected sources 
at z < 1 appear to have much lower accretion rates (e.g. 
iBabic et al.ll2007l ). At decreasing flux levels surveys will be- 
come sensitive to both the very faint sources accreting at low 
redshift and sources accreting at higher rates at high red- 
shift, and therefore the source distribution should become 
double peaked. 

The predicted high redshift peak for the rapid fading 
model is, as expected, more prominent for the scenario where 
more abundant, smaller mass BHs may power quasars. For 
the case with Lpeak,min = 10^ L© almost a quarter of sources 
are predicted to be above z ~ 6 at the lowest flux level 
plotted. The fraction of sources expected above z = 4, 6, 8 
for each of the models at the three flux limits shown in 
Figure 5 is tabulated in Table 1. 



5.2 The integrated X-ray background 

The Cosmic X-ray Background (CXRB) and its unresolved 
component provides an important consistency check for 
models of the number of faint sources. 

We may calculate the contribution of quasars in a red- 
shift band [z, z + dz) to the X-ray background in an observed 
band X by integrating over the QLF, 



dz 



■dz 



dn 



:=/mi„47rD2(2) 



dlogL} 



,dVc 



(14) 

where Dl is the luminosity distance and the limits of inte- 
gration are determined from the limiting sensitivity fmin- 

In Figures 7 and 8 we compare the total CXRB flux 
predicted for the slow fading model, and the rapid fading 
models with low and high Lpeak.min, to the measurement 



of the total CXRB from iMoretti et al.l l|2003l 'l. This contri- 
bution is plotted as a function of the minimum detectable 
flux for sources with redshifts above z — 0.01, 2, 4, , 6, 8 and 
10. We note that the resolved fraction is sensitive to uncer- 
tainties in the absolute value of the CXRB. We restrict our 
analysis of the resolved fraction to energies below 10 keV, 
where variations in the normalis ation between expe riments 
differ by ~ 10 per cent (see, e.g. iMoretti et al.ll2003l ). 

In the 2—10 keV band, we flnd that all models are con- 
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Table 1. Table of the fraction of sources with observed 0.5 — 2 keV 
fluxes above fi^in a-nd redshift above z^i^ for the 3 merger-driven 
models considered. 



sistent with the total measured CXRB. Our models predict 
that ~ 70 - 78 per cent of the total 2-10 keV CXRB is 
due to AGN with fluxes above the best current sensitivity 
level in this band (~ 1.4 x 10~^^ cgs). This result is consis- 
tent with the ~ 80 per c ent resolved fraction measur ed by 
IWorslev et aL I ([2005') and lHickox fc MarkevitchI (12006'). Our 
models remain consistent with the total CXRB at lower flux 
limits. At a sensitivity of 3 x 10~^* cgs - a sensitivity within 
reach of next generation instruments (see Section 5.3) - the 
resolved fraction increases to ~ 87 per cent depending on 
the model. Our models thus still leave room for the very 
hard spectru m sources need to ma ke up the CXRB at ener- 
— dlogLx S^*^^ ^ l|Worslev et al.llioosi ') and/or the population of 
star-forming galaxies expected to contribute significantly at 
low flux levels (e.g. Bauer et aL 2004). We note that our cal- 
culated resolved fractions also depend on the analytic fits we 
have used below z = 2 . Usin g the more recent LDDE model 
from [Silverman et al.l (|2007l ) we obtain lower resolved frac- 
tions (~ 62 — 75 per cent due to sources above the current 
detection level), consistent with the more stringent optical 
selection criteria for this sample. 

In the soft X-ray band (0.5 — 2 keV) , the fraction of the 
total CXRB flux due to sources above the current sensitivity 
level (~ 2.5 x IQ-i'^ cgs) in our mod els is again consi stent 
with the 80 - 90 per cent fou nd bv .Worslev et al.l (|200i) and 



iHickox fc MarkevitchI (|2006l ). However, at fainter flux levels 
the slow fading model overpredicts the total CXRB, again 
reflecting the excess of faint sources at z ~ 2 compared to the 
observed QLF. This contribution again depends on the faint 



Detecting high redshift quasars in X-rays 13 







1 




1 










blow Tading — 






















T 

tn 










en 










cn 












CN 


















> 








\ 


aj 








\ 


o 










T 








\ \ 


CN 


















\ 


O 








\ \ - 










cn 










O 




1 N 


\ 

'J 


\ \ 



-18 -16 -14 -12 
loSio [fmin / erg s"^ cm-^] 



cn — — — 



cn 
O 




Rapid fading- 



cn 
O 



910 [fmin / erg s cm ] 




910 [fmin / erg s cm ] 



Figure 7. Integrated hard (2 — 10 keV) X-ray flux due to AGN for th e mod els in Figure 5. The hatched region shows the estimate of the 
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end behaviour of the fits to the data that we have adopted 
below 2; ~ 2. As mentioned, for the soft X-ray band we have 
chosen to use the LDDE2 model from iMivaii et al] (120011 ) 
which is constructed to reproduce ~ 90 of the total soft 
X-ray background when integrated out to 2 ~ 5. Alternate 
faint e nd extrapolation s may somewhat ease this excess. Us- 
ing the lHasinger et al.l ( 2005) LDDE fit to the LF of type-1 
AGN below z = 2 we find that the slow fading model satu- 
rates, but d oes not overpred i ct, th e soft X-ray background. 
However the lHasinger et al.l 1120051 ) LDDE fit to the soft X- 
ray type-1 QLF only accounts for ~ 35 per cent of the soft 
band CXRB, reflecting the f act th at the sample of type-1 
AGN used in lHasinger et all l|2005h account for only around 
30 per cent of sources at faint and bright flux levels. 



The resolved fractions are naturally lower for the slow 
and rapid fading models with Lpoak.min = lO'^^ Lq (chosen 
to reproduce the soft X-ray QLF at 2 ~ 2), and with this 
choice the slow fading model remains consistent with the 
total CXRB. 

Recently, (jWorslev et"aLl l2005h and 

iHickox fc MarkevitchI 1 20061 ) have pushed the limit for 
the unresolved background further by taking into account 
the stacked emission from galaxies detected with HST and 
IRAC. Indeed we find that if we integrate the flux due 
to sources below the current detection level and compare 
this to the measurement of the unreso lved component as 
derived by iHickox fc MarkevitchI (|2006l ) , the soft and hard 
X-ray components are both too large for the slow fading 
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Figure 10. Redshift distribution of counts in the 0.5 — 2 keV band for the rapid-fading merger driven model and passive evolution models 
as described in text; Merger-driven (solid line), passive evolution with total duty cycle /duty = 0.1 (dashed line) and /duty = 1-0 (dot- 
dashed line) for three flux limits as labelled in cgs units. Note that the lowest flux limit is the anticipated sensitivity for the Generation-X 
mission, the middle for a 1 Ms observation with XEUS, and the right for a point-source within reach of current instruments. The 
horizontal line corresponds to a density of one source per XEUS WFI FOV (7 arcmin)^, which is similar in area to the proposed FOV 
for Generation-X. 



model. Models like this for which sources with relatively 
flat spectra recover the entire unresolved CXRB in the 
soft and hard bands are likely to be in conflict with the 
overall m easured sha pe of the CX RB above 8 keV (e.g . 
IWorslev e t al. 2005; C omastri et al.| [i995: Gil h et al.ll2007l ). 
This consistency check argues further against the slow 
fading model as we have applied it, but we note that a 
flattening of the fading law slope for low mass black holes 
would alleviate the inconsistencies with the CXRB and 
the soft X-ray source counts. In particular the slow fading 
model with I/pcak,min = 10^^ L0 saturates, but does not 
over-predict, the unresolved component of the soft CXRB. 

In a recent paper, ISalvaterra et al.l ()2007l ) predict the 
contribution of high redshift AGN to the unresolved CXRB 
for a merger tree based model for the Eddington limited 
growth of BH seeds. They find that ~ 5 per cent of the 
unresolved 2—10 keV CXRB will be due to sources at z > 6. 
This is similar to the ~ 10 per cent that we predict for 



the rapid fading model with Mbh,min = 3 x 10* Mq, but 
significantly larger than we would predict for models with 
larger Mbh ,min. 

In the next section we focus on the predictions for the 
rapid fading models, since these are consistent with all avail- 
able data, but have included the predictions for the slow 
fading model for reference. 



5.3 Future hard X-ray surveys and the early 
growth of supermassive black holes 

The sensitivities of the deepest X-ray surveys to date are 
photon limited and therefore independent of the background 
intensity. This may not be the case for the effective collecting 
area(s) and angular resolution(s) anticipated for XEUS and 
Constellation-X. Estimates of the point-source sensitivity for 
future instruments therefore potentially become dependent 
on the assumed unresolved extragalactic component (due to 
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quasars and star-forming galaxies) and its flux distribution, 
as well as the effective collecting area and photon extraction 
radius (or resolution). Due to uncertainties in the decompo- 
sition of the X-ray background into contributions from the 
galaxy, star-forming galaxies and quasars, it is somewhat un- 
certain what the confusion limit of these telescopes will be. 
The projected point-source sensitivities for the next genera- 
tion telescope peak for soft X-ray energies. At these energies 
the soft thermal galactic com ponent is a significa nt portion 
of the background (see, e.g. IParmar et al.l 0*9991 ) . however 
the contribution of unresolved point sources may also play 
an important role if the ambitious goals for the resolution 
are not achieved. Since our models predict the contribution 
of quasars to this background, the expected sensitivity of an 
instrument like XEUS to high redshift sources may depend 
on the modelling of faint sources. 

iHasinger et al.l (|2006l ) discuss the anticipated point 
source sensitivity of Constellation-X and XEUS. Assum- 
ing that the observations are confusion limited when there 
are fewer than 40 "beams" per source, combined with es- 
timates of the background due to unresolve d extragalactic 
source s, galactic emission and cosmic-rays, IHasinger et al.l 
(|2006l ) project that a 1 Ms observation with XEUS will yield 
a point source sensitivity of 3 x 10^^* erg s~^ cm~'^, approx- 
imately 200 times better than that of XMM-Newton. This 
estimate presumes that XEUS's goal resolution of 2 arc- 
seconds will be achieved. Achieving this resolution is ex- 
pected to be extremely challenging, and the required resolu- 
tion o f 5 arcseconds is perhaps more realistic (|Hasinger et al.l 
I2OO6I ). For this resolution the sensitivity degrades to ~ 
2 X 10"^'^ erg s~ ^ cm~^ for a 1 Ms int egration. Similarly, for 
Constellation-X fHasinger et al.l (|2006f ) project a point source 
sensitivity of 2 x 10~^^ erg s~^ cm~^ for a 1 Ms observation 
for the Constellation-X goal resolution of 5 arcseconds. 

Assuming a circular beam and applying the above con- 
fusion criterion suggests that XEUS will be able to iden- 
tify point sources with densities below ~ 8.6 x 10* deg~^ if 
it reaches its goal resolution. Comparison with the source 
densities in our models suggests that XEUS should not 
suffer from confusion due to faint AGN at this flux limit 
for the rapid fading models, but would do so below ~ 
10~^^ erg s~^ cm~^ for the slow fading model. The reduced 
sensitivity is not merely attributable to the excess in the 
predicted number of faint sources at 2: ~ 2, since the density 
of sources above z ~ 4 alone are sufficient to stop XEUS 
from reaching its projected sensitivity in this model. This 
suggests that the confusion limit for XEUS may be some- 
what sensitive to the way quasars fade at 2 < z < 6. How- 
ever, we note that the same slow fading model with a higher 
minimum peak luminosity I/pcak,min ~ 10^^ L© would not 
predict a degradation of the estimated point source sensi- 
tivity due to source confusion by faint AGN. Obviously the 
situation would be less favourable if the resolution goal of 
2 arcseconds HEW could not be achieved. If the resolution 
is reduced to the required value of 5 arcminutes, we find the 
density of sources above ~ 5 x 10^^ erg s~^ cm~'^ exceeds the 
density defined by the confusion limit (~ 1.4 x 10* deg~^) 
for the rapid fading model with Lpcak.min ~ 10^ L0, and 
even the rather flat low-flux source density in the rapid fad- 
ing model with I/pcak.min ~ 10^^ would limit the sensitivity 
to sources with fluxes above above ~ 3 x 10^^ erg s~^ cm~^. 



Similar sensitivities would be applicable to Constellation-X 
if it reaches its goal resolution. 

The planned WPfl is designed to have a fleld of view 
(FOV) of (7 arc minutes ) ^ [com pared to (5 arcminutes)^ for 
Constellation-X, lGarciall|2007l )]. An unclustered distribution 
of sources with density ~ 6.7 x 10^ deg~^ will then contain 
one source per WFI FOV. We highlight this density as a 
horizontal line in Figures 9 and 10 for reference. For the 
rapid-fading models, for which observations should not be 
confusion limited down to the goal sensitivity of XEUS, we 
predict significant numbers of sources (more than one per 
WFI FOV) out to 2 ~ 10. To put these sensitivities in con- 
text for models of SMBH growth, in Figure 11 we plot the 
0.5 — 2 keV flux in the observed fra me for BHs ac c reting 
at their Eddington limit (using the iHopkins et al.l (|2007f ) 
2-independent SED) as a function of the source redshift. 
The lowest mass BH that may be seen accreting at its Ed- 
dington limit at z ~ 6 — 10 for a 1 Ms observation with 
XEUS is in the range ~ 10^ - 5 x 10^ Mq. For a model in 
which only BHs larger than ~ 3 x lO'^ M© are active {i.e. 
ipcak.min = lO'^^ Lsun), the rcdshlft limit up to which we 
predict more than one source per XEUS fleld is reduced to 
8 (7) for the slow (rapid) fading model. 

Contrastingly a telescope with 5 arcsecond resolution, 
which as mentioned above may be confusion limited below 
~ 3 X 10~^^ erg s~^ cm~'^, may be sensitive enough to de- 
tect to more than one source per XEUS (Constellation-X) 
FOV up to 2 ~ 7 (2 ~ 6), but detection of signiflcant num- 
bers of sources at higher redshifts would require surveying 
many XEUS/Constellation-X flelds (see Figure 9). This re- 
sult changes little for the case Lpcak.min = lO^'^ Lq, with 
approximately one source per XEUS field out to z ~ 6 pre- 
dicted independently of the fading law. 

The confusion limits we have discussed represent 
the best case f or ea ch model given the assumptions in 
IHasinger et al.l (|2006l ) for the contribution from galactic 
emission and star-forming galaxies, and the true confusion 
limits may therefore be worse if these have been underesti- 
mated. 

Improvements in the resolution of hard X-ray telescopes 
will naturally improve the confusion limit for the point 
source sensitivity (where it is limited by discrete sources), 
particularly for the case of rapid fading in which the source 
counts rise slowly with decreasing flux. NASA's mission con- 
cept Generation-X has a goal resolution of 0.1 arcseconds 
and 100 m^ of collecting area yielding a photometric sen- 
sitivity of ~ 2.2 X 10"^° ergs s"^ cm"^ at 0.1 - 10 keVQ- 
This resolution would certainly avoid confusion due to faint 
quasars, even for the slow fading model where there are 
many faint foreground sources confusing the detection of 
sources at 2 > 5. 

With the proposed sensitivity of Generation-X, unob- 
scured BHs more massive than around a few x 1000 Mq 
would be detectable out to redshifts approaching 15 (see Fig- 
ure 11). Detection of the seeds of SMBHs is one of the prime 
science goals of Generation-X. The number of detectable ob- 
jects at 2 ~ 6 — 15 will obviously depend on the masses of 



ftp://ftp.xray.mpe.mpg.de/peopIe/bol/xeus/XEUS_150108.pdf 
^ http:/ /www.psu.odu/dept/csrp/missions_awarded_genx.htm 
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seed BHs and the growth history of BHs at very high red- 
shifts. 

In Figure 10 we compare the predicted number of de- 
tectable sources for the rapid fading merger-driven model 
[solid Hne) to the passive evolution models with /duty = 0.1 
[dotted line) and /duty = 1-0 [dot-dashed line) for three val- 
ues of the point-source sensitivity. The left panel shows the 
evolution in the density of sources with fluxes above the sen- 
sitivity goal of Generation-X. The space density of observ- 
able sources at z ~ 15 depends very sensitively on the model 
assumed. The FOV of Generation-X is expected to be at 
least (5 arcminutes)^, with fiel ds as big as (15 arcmin utes)^ 
considered in the design study (| Windhorst et aLll2006l ). Here 
we assume the FOV of Generation-X is the same as for 
XEUS for ease of comparison. For the merger driven model, 
with a minimum active BH mass of ~ 3 x 10** M© we predict 
significant numbers of sources, still around 10 per FOV, out 
to z = 15. The situation is even more promising if the BHs 
are evolving passively with a small duty cycle, with up to 
~ 100 sources per field up to z = 15. If the BHs are evolving 
rapidly there is still more than one source per field out to 
2 ~ 11, but the density drops very rapidly with increasing 
redshift for a scenario in which all BHs are growing via ac- 
cretion at their Eddington limit. At the goal sensitivity of 
XEUS, the maximum redshift at which there is one source 
predicted per FOV also varies widely; ranging from z ~ 8 
for continuous growth to z ~ 15 for slowly growing BHs. 



6 SUMMARY AND CONCLUSIONS 

We have presented a hybrid model for the redshift evolution 
of the X-ray emission connected with the fuelling of super- 
massive black holes with the aim of assessing the prospects 
of detecting quasars via their X-ray emission at z > 6. At 
z < 2 we have used the observed X-ray luminosity func- 
tions, at 2 < z < 6 we have adopted a CDM merger-driven 
model for the evolution of the emission from supermassive 
black holes combined with a slow and a rapid fading law and 
we have explored a range of assumptions for the growth of 
supermassive black holes at z > 6. 

For appropriate choices for the efficiency of black hole 
formation in a dark matter halo and the characteristic 
quasar lifetime, our model is in good agreement with the 
observed optical and soft and hard X-ray quasar luminos- 
ity functions at < 2; < 6 for both the slow and the rapid 
fading law. The only disagreement that occurs is for the 
faint end of the X-ray luminosity functions where our merger 
model combin e d with the slow fading model suggested by 
iHopk ins et al] (|2005al ) based on detailed numerical simula- 
tions predicts too many faint objects. It consequently also 
overpredicts the soft X-ray background suggesting that the 
slow fading model of [Hopkins et al. ( 2005a) cannot extend 
to black hole masses much below lO'^ Mq. With a more 
rapid fading law our models are well within the limits of 
the integrated CXRB, leaving room for the population of 
star-forming galaxies expected at low fiux densities, and the 
population of Compton-thick sources needed to explain the 
overall spectral shape of the CXRB. 

The main aim of our study is to assess the prospects 
of planned and proposed X-ray missions like Constellation- 
X, XEUS and Generation-X to study the build-up of black 



holes at z > 6. With a point source sensitivity of 3 x 
10~^^ ergs s~^ cm~^ for a 1 Ms exposure, Constellation- 
X will detect significant numbers of black holes at 2: ~ 5 — 6 
but will probably not yet reach the necessary sensitivity to 
detect significant numbers of sources at 2 > 6. This conclu- 
sion holds independently of our assumptions for the growth 
of BHs at z > 6. For XEUS with its anticipated ten times 
superior point source sensitivity, the prospects for detecting 
sources at 2 > 6 is much more favourable. Up to 17 per cent 
of the approximately 100 sources in the 49 square arcminute 
FOV expected for the rapid fading model with a minimum 
active black hole mass of ~ 3 x 10* M© are expected to 
be at 2 > 6. If our merger driven model can be extrapo- 
lated to 2 > 6 or if black holes grew with a duty cycle of 
about 10 per cent, XEUS would detect significant numbers 
of black holes with a rather fiat redshift distribution out to 
2 ~ 10. In either case, observable X-ray emission would have 
to accompany the growth of supermassive black holes start- 
ing from seed black holes of > 10* Mq or smaller. For the 
more remote prospect of Generation-X, with an anticipated 
point source sensitivity yet a factor 100 better again, the flat 
redshift distribution could extent to 2 ~ 15 and beyond. 

If black holes grew much faster at 2 > 6 as in our Ed- 
dington limited growth model with a duty fraction of unity, 
which would be required for Eddington limited growth of 
the most massive black holes from stellar mass seed black 
holes, neither XEUS nor Generation-X would detect many 
black holes at 2 > 10. However, while such rapid growth 
may occur and may indeed be necessary for the most mas- 
sive black holes at 2 ~ 6, it appears unlikely that this will be 
the norm for the majority of black holes. The prospects that 
XEUS (and Generation-X) will unravel the growth history 
of black holes at 2 > 6 and beyond are thus excellent. 
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